arXiv: 1503.0254lv2 [cond-mat.mes-hall] 15Jul2015 


Fingerprints of disorder in graphene 


Pei-Liang Zhao, 1 Shengjun Yuan, 2 EI Mikhail I. Katsnelson, 2 and Hans De Raedt 1 

1 Department of Applied Physics, Zernike Institute for Advanced Materials, 

University of Groningen, Nijenborgh 4, NL-9747AG Groningen, The Netherlands 
2 Radboud University, Institute for Molecules and Materials, 

Heijendaalseweg 135, 6525A J Nijmegen, The Netherlands 
(Dated: July 17, 2015) 

We present a systematic study of the electronic, transport and optical properties of disordered 
graphene including the next-nearest-neighbor hopping. We show that this hopping has a non- 
negligible effect on resonant scattering but is of minor importance for long-range disorder such as 
charged impurities, random potentials or hoppings induced by strain fluctuations. Different types of 
disorders can be recognized by their fingerprints appearing in the profiles of dc conductivity, carrier 
mobility, optical spectroscopy and Landau level spectrum. The minimum conductivity 4e 2 //i found 
in the experiments is dominated by long-range disorder and the value of 4e 2 /7r h is due to resonant 
scatterers only. 

PACS numbers: 72.80.Rj; 73.20.Hb; 73.61.Wp 


I. INTRODUCTION 

The dominant source of disorder which limits the 
transport and optical properties of graphene is still un¬ 
der debate. Different mechanisms have been proposed 
and investigated intensively, including charged impuri¬ 
ties, random strain fluctuations and resonant scatter¬ 
ers (for reviews see Refs. [lp). Early on, charged im¬ 
purities (Cl) have been recognized as the dominate dis¬ 
orders due to graphene’s unusual linear carrier-density- 
dependent conductivity. However, this mechanism does 
not explain the experimental observations that the trans¬ 
port properties of certain samples are not sensitive to 
the substrate screening— On the other hand, strain 
fluctuations (SF) induced e.g. ripples can be alternative 
scattering mechanism—; they can be also responsible for 
charge inhomogeneities, that is, electron-hole puddles&£. 
There is experimental evidence, based on the correlation 
between the carrier mobility and the width of the resis¬ 
tance peak around charge neutrality, that the long-range 
disorder potential (LRDP) due to SF could be the dom¬ 
inant source of disorder in high-quality graphene on a 
substrate^. In addition, the SF modulate the electron 
hopping energies between different atomic sites, induc¬ 
ing the long-range disorder hopping (LRDH), leading 
to the appearance of the (pseudo) vector potential^. 
Another common source of disorder are resonant scat¬ 
terers (RS) such as chemical species like hydrogen or 
organic groups, which also lead to a sublinear carrier- 
density-dependent conductivity and a minimum conduc¬ 
tivity plateau around the neutrality poin t 10 : 11 . 

Besides the transport properties, an important part 
of our knowledge about the electronic properties de¬ 
rives from the optical spectroscopy measurement s 1 : 1 2 . 
Infrared spectroscopy experiments allow for the con¬ 
trol of interband excitations by means of electrical gat- 
in g 13 : 14 . For doped pristine graphene with nonzero chem¬ 
ical potential pLp, the optical conductivity is a step func¬ 
tion a (co) = ctq© (w — 2/ij?) at zero temperature due to 


Pauli’s exclusion principle. However, there are experi¬ 
mentally observed background contributions to the opti¬ 
cal spectroscopy between 0 < u < 2^ f 14 ' 15 , which are 
due to the extra intraband excitations introduced by dis¬ 
order or many-body effect s 14 : 16- —. This opens the pos¬ 
sibility to identify the source of disorder via the optical 
measurements. 

Previous theoretical investigation of disorders are 
mainly based on models without considering the next- 
nearest-neighbor (NNN) hopping t'. The breakdown of 
electron-hole symmetry resulting from f / 0 shifts the 
position of Dirac point from zero to 3t — >— Recent quan¬ 
tum capacitance measurements indicate that the value of 
t' is about 0.3eV— , consistent with the values obtained 
from the density functional calculations. It is generally 
thought that t' has relatively weak effects on the phys¬ 
ical properties of graphene at low energies^ : 22 : 27 : 28 . In 
the present paper, we study the electronic, transport and 
optical properties of graphene with different types of dis¬ 
orders including NNN. We show that if has a negligible 
effect in combination with long-range disorder such as Cl, 
LRDP and LRDH, but changes the physics dramatically 
when RS are present. Different sources of disorder can be 
identified via their fingerprints in the common measur¬ 
able quantities, such as dc conductivity, carrier mobility, 
optical spectroscopy and Landau level spectrum etc. We 
will use these fingerprints to demonstrate the dominant 
disorder source in several well-known experimental mea¬ 
surements. The paper is organized as follows. In Section 
II we gives a description of the tight-binding Hamilto¬ 
nian of single layer graphene including different types of 
disorders and NNN. In section III and IV, we discuss the 
effect of different disorders on the transport and opti¬ 
cal properties of graphene. Then, we study the Landau 
level spectrum and quantum capacitance in the presence 
of perpendicular magnetic field in section V. Finally, a 
brief discussion and conclusion, including a list of domi¬ 
nant disorder sources in several experiments, is given in 
section VI. 
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FIG. 1: (Color online) The dc conductivity as a function of 
carrier density n e for disordered graphene. Left panels show 
the results without the NNN hopping t ', and right panels with 
t' = O.lt. For Cl, we use k = 6 of hexagonal-boron nitride 
as a typical value of dielectric constant for graphene on a 
substrate. The use of other k for different substrate such as 
SiC >2 does not change the results quantitatively. Here 0.01% 
disorder corresponds to a concentration of 3.82 x 10 11 cm -2 . 


II. MODEL AND METHOD 


We consider disordered graphene described by the 
tight-binding (TB) Hamiltonian 

H = - hA 0 ! - H t'i,j C i C j + Vi Ci Ci ' ( 1 ) 

i,j i,j i 





n e (lo'W 2 ) n e (l0 12 cm' 2 ) 



FIG. 2: (Color online) The carrier mobility as a function of 
carrier density n e for disorder graphene with t'. 


where the first sum is taken over nearest neighbors and 
the second one is over next-nearest neighbors. 

For Cl, we consider randomly distributed point-like 
charges at the center of a hexagon of the honeycomb lat¬ 
tice (r fc )Sa, which introduce the Coulomb energy Vi = 
sign (&) e 2 / (k|d — r fc|) at eac h site i, and the screen¬ 
ing effect due to the substrate is taken into account by 
using the dielectric constant k of the substrate. Here, ac¬ 
cording the values of sign (k) we consider three types of 
Cl: (1) CI° for randomly distributed positive or negative 
potential caused by charges that the whole sample holds 
the electric neutrality, (2) CI + for only positive potential 
and (3) CI _ for only negative ones. 

For LRDP, the on-site potential Vi follows a corrected 
Gaussian profile which varies smoothly on the scale of 
lattice constant as v% = Yh k £4 exp[— |r» — rfe| 2 /(2d 2 ))^, 
where r k is the k — th Gaussian centers which are ran¬ 
domly distributed on the lattice with probability p v , U k 
represents the strength of the local potential and is uni¬ 
formly random in the range [—A„,A„], and d is inter¬ 
preted as the effective radius. We use A v = t and d = 5a 
to represent the long-range Gaussian potential. Here 
a « 1.42A is the carbon-carbon distance in the single¬ 
layer graphene. 

The LRDH is introduced in a similar way as LRDP 
except that the nearest-neighbor hopping parameters 
are modified according a correlated Gaussian form via 
Uj = t + J2k T k ex P[— |D + r,• - 2r fc | 2 / (8df)], where T k , 
dt and pt have similar meanings as in LRDP, and we 
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choose At = 0.25 1 and dt = 5a— We want to emphasis 
that, although the amplitude (A) and radius (d) of the 
Gaussian profile in the LRDH and LRDP are free pa¬ 
rameters that can be turned in the tight-binding model, 
the numerical results show little quantitative difference as 
long as these parameters are of the same order as the cho¬ 
sen values. In general, an increase (decrease) of the am¬ 
plitude or radius is equivalent to an increase (decrease) 
of the disorder concentration. 

The hydrogen-like RS is described by the Hamiltonian 

Hrs = (dj c i + h.c) n i 30 i 31 , where V is the hop¬ 

ping between carbon and adatom. We consider the lim¬ 
iting case with V —> oo, i.e., the electron at the impurity 
site is completely localized such that the resonant scat- 
terer behaves like vacancy— In our calculations, we use 
t ~ 2.7 eV and t' = t /10 for the nearest and next-nearest 
neighbor hopping parameters, respectively. The spin de¬ 
gree of freedom contributes only through a degeneracy 
factor and, for simplicity, has been omitted in Eq. ©■ 

The calculations of the electronic and optical prop¬ 
erties are performed by the tight-binding propagation 
method (TBPM) 11 ’ 31- — , which is based on the numer¬ 
ical solution of the time-dependent Schrodinger equation 
and Kubo’s formula. The advantage of this method is 
that all the calculated quantities are extracted from the 
real-space wave propagation without any knowledge of 
the energy eigenstates. Furthermore one can introduce 
different kinds of (random) disorder by constructing the 
corresponding TB model for a sample scaling up to mi¬ 
crometers. For more details about the numerical methods 
we refer to Refs [2611311 . The simulated graphene sample 
contains up to 8192 x 8192 atoms subject to periodic 
boundary conditions. 


III. TRANSPORT PROPERTIES 

We first consider the carrier-density-dependence of the 
microscopic conductivity a (n e ) for disordered graphene. 
The microscopic (or semi-classic) conductivity is calcu¬ 
lated from the diffusive region of the charge transport, 
i.e., when the time-dependent diffusion coefficient reaches 
its the maximum^ 3 -, and it is comparable to the con¬ 
ductivity extracted from the field-effect measurements. 
In TBPM, the microscopic conductivity at an energy E 
is calculated by using the Kubo formul a 31 ’ 33 

a (E) — max^P [ dt Re \e~ lEt {ip\ Je im J |E)1 , 

T “ Jo 

( 2 ) 

where |<p) is a normalized random state, | E) is the nor¬ 
malized quasi-eigenstate^i, J is the current operator, is 
the sample area, and p ( E) is the density of states (DOS) 
calculated vi a 31 ’ 32 

1 r°° 

P(E) = — e im {v\v{t))dt. (3) 


The measured field-effect carrier mobility is related to 
the microscopic conductivity as p (E) = a ( E ) /en e (E), 
where the carrier density n e is obtained from the integral 
of DOS via n e (E) = p (e) de. 

From the results shown in Fig. [I] we see that (1) in¬ 
cluding t' has negligible effects for Cl, LRDP and LRDH, 
but the results for RS change dramatically. In the pres¬ 
ence of RS, there is a strong electron-hole asymmetry in 
the carrier-density-dependence of dc conductivity. This 
is due to the fact that the impurity band created by RS 
is shifted from the Dirac point to the hole side2I, in¬ 
troducing strong electron-hole asymmetry at low ener¬ 
gies; (2) as a consequence of this shift the conductiv¬ 
ity plateau around the neutrality point is also shifted 
to the hole side, with an impurity-concentration depen¬ 
dent height and width (for very small concentration of 
RS, there is just as a kink instead of a plateau, see the 
point indicated by an arrow in Fig. [ijh)) i These fea¬ 
tures can be observed in graphene if the concentration 
of generic RS is increased by exposing the material to 
atomic hydrogen—. (3) a ( n e ) exhibits a sublinear depen¬ 
dence for small concentration for all types of disorders, 
except for the hole side in the presence of RS; (4) For 
LRDH, a (n e ) is insensitive to the changes of the disorder 
concentration ( p t ); (5) No matter whether t' is nonzero 
or not, linear-dependent cr (n e ) appears only in Cl with 
large concentration of indicating that Cl is the 

dominant source of disorder in the experimental samples 
which show clearly the linear carrier-density-dependent 
conductivity (such as K151 in Ref.—, and Potassium 
doped samples in Ref—, etc.), agree with the theoret¬ 
ical prediction that a (n e ) oc n e ; (6) The electron-hole 
asymmetry appears also for larger concentration of Cl 
if there is only one types of charge resource (CI + and 
Cl - ). However, this asymmetry is different from the one 
due to RS in two aspects: first, for Cl there is no kink 
or plateau in the profile; second, the conductivities on 
both electron and hole sides decrease significantly with 
larger concentration of Cl; (7) Only in the case of CI + the 
conductivity on the electron side is smaller than on the 
hole side with the same concentration of carrier density, 
which is a unique signature of CI + . This is in concert 
with experiment result s 15 ’ 39 . 

The field-effect carrier mobility p can be calculated 
from the conductivity and carrier-density through p = 
a / en e . In the following we show only the results with 
non-zero t'. From the results presented in Fig. [2j we see 
that (1) the carrier-dependence of mobility p ( n e ) is very 
similar for CI° and LRDP; (2) for LRDH, p (n e ) is insen¬ 
sitive to the disorder strength; (3) electron-hole asymme¬ 
try appears for CI + , Cl - and RS, but only in the case of 
CI + the electron mobility is smaller than the hole for the 
same concentration of carrier density; (4) for RS, the mo¬ 
bility on the electron side is insensitive to the impurity 
concentration, and its value can be one order of magni¬ 
tude larger than the value on the hole side. For example, 
considering a RS concentration of n x = 0.025%, the elec¬ 
tron mobility at carrier density 5 x 10 12 cm -2 is about 
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FIG. 3: (Color online) The optical conductivity as a function 
of energy for disordered graphene with /j,f = 0.1 1 and t' = 
O.lt. Here ao = ne 2 /(2h) is the universal optical conductivity 
of graphene. All along the work the temperature of optical 
calculation is T = 45A', the same as in the experiment of 
Ref[Il. 

3,000 (cm 2 V~ 1 s~ 1 ) but the hole mobility for the same 
carrier density is only ~ 300. This significant one order 
difference of the electron and hole mobility is a unique 
signature of RS; (5) with RS present, on the hole side, the 
carrier-density-dependent mobility is not monotonic and 
/i (n e ) reaches a minimum at the density corresponding 
to the tail of the conductivity plateau. However with RS 
present and t' = 0, the drop of mobility at the minimum 
is one order of magnitude larger than the experimental 
result. 

The minimum conductivity cr m ; n at the Dirac point is 
of the order of 4 e 2 /h for all types of long-range disorders 
with t! = f/10. The values of er m in hr Cl and LRDP do 
not depend on t ', but change with the disorder strength 
such that larger concentration of disorder leads to larger 
values of cr m i n - This is due to the fact that the increase 
of potential sources in Cl and LRDP will increase the 
DOS at the fip, leading to more states which can con¬ 
tribute to the transport. This may also explain the exper¬ 
imental observations in Ref. (39| and Ref. (4lj in which the 
low mobility does not necessary correspond to a smaller 
value of (J m i n . For LRDH, the value of cr min for t' = 0 
is about two time larger than the value for t' = t/10, 
but both are insensitive to the disorder strength. For RS 
and t' = 0, (7 m i n is of the order of 4e 2 /7r h, independent on 
the impurity concentration , but if t' = t/10, er m in 


from being of the order of 4 e 2 /h at small n x to Ae 2 /ph 
when n x > 0.1%, consistent with the numerical results 
of Ref. M (data not shown). Thus we conclude that 
our results indicate that the minimum conductivity 4e 2 / h 
found in the experiments is dominated by long-range dis¬ 
order but that the value of Ae 2 /p h is due to RS only. It 
is worth to mention that our consideration does not take 
into account the effects of weak (anti) localization which 
can change the behavior of conductance at very large 
distances^, due to energy smearing in our calculations. 
The latter works as dephasing. At the same time, this 
dephasing can be physical for real samples. 

IV. OPTICAL SPECTROSCOPY 

The optical conductivity is calculated by using the 
Kubo formula^ within TBPM4I as (omitting the Drude 
contribution at u = 0) 

e -/3^ _ i r°° 

a (ui) = lim -—— / e _£t sin uit 

e- J-0+ Wl2 J q 

x2 Im (tp\f (H) J (t) [1 - / (H)] J\ip) dt , 

(4) 

where /3 = 1 /ksT is the inverse temperature, f (TL) = 
1/ + l] is the Fermi-Dirac distribution oper¬ 

ator. Similar as for the transport properties, our nu¬ 
merical calculations show that t' has negligible effects 
on the optical properties of disordered graphene, ex¬ 
cept if RS are present. In general, disorder introduces 
new states which could contribute to the extra intra¬ 
band excitation o 14 ! 16- — , and therefore enhances the op¬ 
tical conductivity below 2fj,p , which might explain the 
observed background contribution in the optical spec¬ 
trum for 0 < ui < 2/tj i 14 i 15 . This is confirmed by the 
optical conductivity of disordered graphene calculations 
shown in Fig. [3] For disordered graphene with Cl (in¬ 
cluding CI°, CI + and CI _ ) there is a strong enhancement 
of the optical conductivity below 2 /ip and the enhanced 
spectrum forms a plateau with disorder-dependent min¬ 
imum conductivity. For LRDP, there is in addition a 
disorder-dependent plateau in the optical spectrum be¬ 
low 2fip, which is much wider that the one due to CI. 
For LRDH, the enhancement of the optical conductivity 
is much smaller than for other types of disorders. For 
RS and t' = 0, a disorder-dependent peak appears at 
ui ss /ip, which is due to the enhanced excitations of the 
midgap states at the Dirac point. This peak disappears 
for t' = t/ 10, and instead, a disorder-dependent narrow 
plateau appears. 

In practice, instead of varying the disorder concen¬ 
tration, it is easier to change the chemical potential by 
applying an electrical potential to a gate. In order to 
compare to the experimental data of the spectroscopy 
measurements 14,15 quantitatively, we plot in Fig. [I] the 
best fit of the optical conductivity for different chemical 
potentials ranging from 0.05f to O.lt (since the results 
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FIG. 4: (Color online) The optical conductivity as a function of energy for disordered graphene (a) LRDP, (b) Cl, and (c) 
RS with t'. The disorder concentrations are determined via the best fit to the experimental results of optical spectroscopy— 
The chemical potential fiF changes from 0.05f to O.lt. The inner panels show the corresponding carrier mobility for the same 
concentration of disorder. The dashed line in (c) is a guide to eye separating two region in which the spectrum changes 
differently in the presence of RS. 


of CI°, CI+ and Cl" are similar, we present here only 
the case of CI°). The disorder concentrations shown in 
Fig. 0] are determined by matching the minimum value of 
the optical conductivity plateau to the one observe d 14 : 15 , 
yielding a p iateau of the order of O.luo for p,p ~ f/10. The 
best match of the disorder concentrations from our simu¬ 
lations is p v = 0.01% for LRDP, nc = 0.025% for Cl and 
n x = 0.025% for RS. A direct comparison of the profile 
of the spectrum between our simulations and the experi¬ 
ments in Ref. [mi indicates that LRDP fits best to the 
experiments. In Ref. Il4l . the carrier mobility measured 
for the same device is as high as 8, 700cto 2 U _1 s _1 at car¬ 
rier densities of 2 x 10 12 cro -2 , and the LRDP also gives 
the highest mobility that it can reach ~ 3,000. For Cl, 
/i ~ 1500, and for RS the mobility is even smaller: for 
electrons It is ~ 1,000 and for holes ~ 300. Therefore we 
conclude that the background contribution of the optical 
conductivity below 2//jr as observed in Ref. M should be 
due mainly to the presence of LRDP. 

V. LANDAU LEVEL SPECTRUM 

Finally we consider the electronic properties of 
graphene under a perpendicular magnetic field (B = 
50T). The Landau quantization of the energy levels leads 
to separated peaks, as shown in Fig. [5] In the presence of 
disorder, the peak amplitudes of the Landau levels (LL) 
are reduced and the peaks become broader, except for 
LRDH in which the influence of disorder is much weaker 
than for other types of disorders. The peak profiles de¬ 
pend on the different sources of disorder. In general, for 
long range disorder, the peak is still symmetric along its 
center, but for RS, the changes are mainly restricted on 
the side with higher energy. Furthermore, the LL spec¬ 
trum exhibits electron-hole symmetry for CI° and LRDP, 
but becomes asymmetric for CI + , CI _ and RS. Espe¬ 
cially, there are two small peaks around the first Landau 


level on the hole side shown in Fig. [5jd) , which has the 
same origin as for the zero LL peaks, induced by RS^I. 
The differences that appear in the LL spectrum also ap¬ 
pear in quantum capacitance measurement, as the inverse 
of the latter is proportional to DOS—Therefore, we 
also expect a huge effect of RS on the asymmetric quan¬ 
tum Hall conductivity, a topic for future research. 

The quantum capacitance C qi which is defined as 
C q = pe, can be extracted experimentally from the to¬ 
tal capacitance C and the geometry capacitance C g via 
1 /C q = 1/C — 1 /C g . In Fig. [6] we show the carrier de¬ 
pendence of 1/p, which is proportional to 1 /C q , for dif¬ 
ferent types of disorders under the same magnetic field 
(B = 50T). Due to the presence of disorder, the peak 
amplitudes decrease significantly except for the LRDH, 
in which the influence of random hopping is negligible. 
The change of the spectrum profile for each type of disor¬ 
der has similar feature deduced from the corresponding 
DOS. Furthermore, some characters become even more 
clear in the spectrum of 1/p. For example, the electron- 
hole asymmetry appeared in the presence of single-type 
charge impurities (CI + or CI _ ) is very special: the slopes 
of the peaks on the hole and electron sides point to the 
same direction, depending on the sign of Cl (see a zoom 
of the first two peaks in Fig. This unique feature has 
also been observed in the experiments—1 

VI. DISCUSSION AND CONCLUSION 

We have studied the effects of different types of disor¬ 
ders on the electronic, transport and optical properties 
of graphene. By comparing the results with and without 
the NNN hopping, we find that the NNN hopping has 
negligible effect in combination with long-range disorder 
such as Cl, LRDP and LRDH, but that it changes the 
physical properties dramatically if RS are present. In the 
latter case, we find that 1) there is an extra conductiv- 
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FIG. 5: (Color online) Density of states as a function of energy for disordered graphene in the presence of a uniform 

perpendicular magnetic field (B = 50T). 


ity plateau on the hole side, with a value larger than 
the minimum conductivity at the neutrality point; 2) 
the carrier-density-dependent mobility does not always 
drop with larger carrier density but instead, it reaches 
a minimum at the edge of the conductivity plateau. 3) 
a strong electron-hole asymmetry appears in the carrier- 
density-dependent transport properties and Landau level 
spectrum; 4) the minimum conductivity at the shifted 
Dirac point is no longer a constant, but drops to 4e 2 /7r h 
when the impurity concentration is larger than 0.05%. 
For long-range disorder, the minimum conductivity for 
t' = f/10 is of the order of 4e 2 /h and increases with larger 
disorder concentration for Cl and LRDP, but remains the 
same for LRDP. The mobility always becomes smaller 
with larger concentration of disorder, however, the min¬ 
imum conductivity does not follow the same rule, con¬ 
sistent with the transport measuremen t 39 ! 41 . For doped 
graphene, the presence of disorder introduces extra exci¬ 
tations below 2 up but the profiles of the optical spectra 
are different for different types of disorders. 


As an example of using the fingerprints discussed in 
the main text, we collect the dominant source of disorder 
in several well-known experiments and list them in Table 
I. Different types of disorders such as Cl (including CI°, 
CI + and CI _ ), LRDP and RS have been identified in 
different experiments, except for the LRDH which has 
been proved to have negligible influence to the electronic 
properties. The results obtained in Table I also suggest 
the dominant source of disorder may vary from sample 
to sample. 


In summary, we suggest that the different but charac¬ 
teristic features that appear in the calculated electronic, 
transport and optical properties can be used as finger¬ 
prints to identify the dominant sources of disorder in 
graphene. 














































































































































FIG. 6: (Color online) The reciprocal of DOS as a function of carry density n e for disordered graphene in the presence of a 
uniform perpendicular magnetic field (B = 50T). 
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Ref. 

disorder 

Fingerprints 


LRDP 

Symmetrical a (n e ) in Fig.2 (a), the minimum conductivity plateau in Fig. 2 (b), and the relation of 

mobility /x versus n* in Fig. 2 (c). 


id 

RS 

Asymmetrical a (n e ) of the blue and red curves in Fig. 2 (a). 


3 

LRDP 

A plateau in the doped optical spectroscopy in Fig. 2 (b), together with the corresponding relatively 

high mobility. 

Ua 

CI+ 

A narrow plateau in the doped optical spectroscopy, together with a shift of the minimum 
conductivity to the electron side in Fig. 1. 

Hr 

CI+ 

The electron mobility is smaller than the hole one in Fig. 2 for samples K130, K145, K151; The 
minimum conductivity shifts to the electron side in Fig. 3. 

mr 

Cl" 

The hole conductivity is smaller than electron one and the minimum conductivity shifts to the hole 

side in Fig. 2. 


m 

CI+ 

The hole conductivity is larger than electron one in Fig. 1 (the sample before annealing). 


49 

Cl 

A narrow plateau in the doped optical spectroscopy in Fig. 3 (b). 


TABLE I: List of the dominant disorder source in different experimental samples, identified by using the fingerprints appeared 
in the transport or optical properties. The figures indicated in the table are these in the corresponding reference. 




n e (io 12 cm 2 ) 


FIG. 7: (Color online) The reciprocal of DOS as a function 
of carry density n e for disordered graphene in the presence of 
a uniform perpendicular magnetic field (B = 50T). 
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